device,decomposed=0
u = dblarr(64,64)
loadct,5
openr,1,'fort.9',/f77
readu,1,u

close,1
end
